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Abstract 

We present a measurement of the branching fraction for the decay af (1260)7r^ with 

af (1260) Tr^TT^TT^ using a data sample containing 535 x 10^ BE pairs collected with the Belle 
detector at the KEKB asymmetric-cncrgy e~^e~ collider operating at the T(4S') resonance. We 
measure the branching fraction B{B^ af (1260) iT^)B{af {1260) -> vr^vr^vr^F) = (14.9 it 1.6 ± 
2.3) X 10~^, where the first and second errors are statistical and systematic, respectively. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 
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The decay B'^ — > af (1260)7r''' fl*] proceeds through h ^ u transitions, hence its time- 
dependent CP violation is sensitive to the Cabibbo-Kobayashi-Maskawa (CKM) [2| angle 
02 01 • The first attempt to search for this decay mode was made by the CLEO collabo- 
ration; an upper limit of 490 x 10~^ was obtained at the 90% confidence level (C.L.) 
Recently, the BaBar collaboration reported the first measurement of the branching fraction 
of — s> af (1260)7r^ assuming a 50% branching fraction for af (1260) — > tt^tt^tt^ decay: 
(33.2 ± 3.8 ± 3.0) X 10"^, where the first (second) error is statistical (systematic) error p]. 
The measured branching fraction is compatible with the prediction in Ref. j^. BaBar also 
measured the time-dependent CP- violating parameters [7|. The CP- violating parameters 
provide an effective 02 value that can be shifted from the CKM angle 02 due to the contri- 
bution from b —>■ d "penguin" decay amplitudes. The value of 02 can be extracted using the 
method proposed in Ref. 

In this report, we present a measurement of the branching fraction of the decay P° 
af (1260)7r^ with af (1260) — > vr^vr^vr^ based on a data sample containing N^-^ = (535 ± 
7) X 10® BB pairs. It is known that the af (1260) decays into p'^n^ via both S and D-waves, 
and cr(600)7r^ [9|]. In this analysis, we use a Monte Carlo (MC) simulation with S-wave 
af (1260) p'^TT^ decays to estimate the nominal detection efficiency. We then assign 
a systematic error on the detection efficiency from the contributions of af (1260) decays 
into D-wave p^ir^ and (T(600)7r±. The non-resonant contributions from P° — ^ p°7r+7r and 
7T~^7r~n~^7r~ decays are neglected in the nominal fit. We estimate the possible contribution 
of these modes from a fit to the three pion mass distribution and assign a corresponding 
systematic error. 

The data sample was collected with the Belle detector lOj at the KEKB asymmetric- 
energy e~^e~ (3.5 on 8 GeV) collider [Tf]] operating at the T(4S') resonance. The Belle 
detector is a large-solid-angle magnetic spectrometer that consists of a silicon vertex detector 
(SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold Cherenkov 
counters (ACC), a barrel-like arrangement of time-of- flight scintillation counters (TOF), 
and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located inside a 
superconducting solenoid coil that provides a 1.5 T magnetic field. An iron fiux-return 
located outside of the coil is instrumented to detect mesons and to identify muons (KLM). 
A sample containing 152 x 10® BB pairs was collected with a 2.0 cm radius beampipe and 
a 3-layer silicon vertex detector, while a sample with 383 x 10® BB pairs was collected with 
a 1.5 cm radius beampipe, a 4-layer silicon detector, and a small-cell inner drift chamber. 

We reconstruct P° af (1260)71^^ candidates from combinations of four charged tracks 
originating from the beam interaction region. The tracks are required to be consistent with 
a pion hypothesis based on the particle identification (PID) information from the ACC and 
the dE/dx measurements in the CDC. Positively identified electrons are rejected. 

An af (1260) candidate is reconstructed from three charged pions tt^tt^tt'^ with invariant 
mass in the range from 0.8 GeV/c^ to 1.8 GeV/c^. We require at least one pair of oppositely 
charged pion candidates satisfy the condition 0.55 GeV/c^ < m^+^- < 1.15 GeV/c^, where 
m^+^- is the invariant mass of the pair. 

To construct B'^ candidates we combine af (1260) candidates with a bachelor pion having 
momentum in the range 2.2 GeV/c < pbach < 2.7 GeV/c in the T(45') center-of-mass system 
(CMS). The P° candidates are identified using the energy difference /S.E = — -Ebeam ^^nd 



the beam energy constrained mass Mbc = y (-^beam)^ ~" where -Ebeam is the CMS beam- 

energy, and E'^ and p*^ are the CMS energy and momentum of the P° candidate. We select 
the P° candidates in the region 5.20 GeV/c^ < Mbc < 5.30 GeV/c^ and \AE\ < 0.12 GeV. 
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We explicitly eliminate three charm B decay modes that peak in the Mbc signal region: 
5° ^ D-71+ with D' -> K+n~7i', S° ^ D*-7i+ with D*- Wn- and W K+ti-, 
and -B° J/ipK*^ with J/ip ^ fi^fi^ and ^ f^+vr". We exclude candidates in the 
mass ranges: 1.86 GeV/c^ < M^- < 1.88 GeV/c^, 2.00 GeV/c^ < M^*- < 2.05 GeV/c^, 
and 3.05 GeV/c^ < Mj/^ < 3.15 GeV/c^. 

The dominant background for this analysis comes from continuum events of the type 
e~^e~ qq {q = u,d,s,c). To suppress this background, we construct likelihood functions, 
Ls and Lbg, for the signal af (1260)7r^ decay and the continuum background, re- 

spectively. The likelihood functions are formed as a product of probability density functions 
(PDF) of a Fisher discriminant formed from event shape variables [13] and the angle of the 
B candidate flight direction in the CMS with respect to the beam axis. We also make use 
of the variable r provided by a flavor tagging algorithm that identifies the flavor of the 
accompanying B^ in the T(4S') — > B^B^ decay. The parameter r ranges from r = for no 
flavor discrimination to r = 1 for unambiguous flavor assignment. The data sample is di- 
vided into six r intervals. Since the separation of the continuum background from the signal 
depends on r, we determine the requirement on the likelihood ratio 71 = Ls/ {Ls + Leg) for 
each r bin by optimizing the expected sensitivity using signal MC events and events in the 
sideband region Mbc < 5.26 GeV/c^. 

We find that on average 2.9 B^ — > a]'^(1260)7r'^ decay candidates are included for each 
event selected from the data. The B^ candidate having the largest TZ value is chosen in 
the case of multiple candidates in an event. After the best candidate selection, we find 
that 19.9% of the signal events in the MC simulation are incorrectly reconstructed self-cross 
feed (SCF) events; for those events at least one charged pion track is replaced with one 
from the accompanying B meson decay. In 99.6% of the SCF events one or more tracks 
from the a^(1260) — >■ TT^TT^TT^ has been replaced; in the remaining 0.4% a bachelor pion 
was replaced. For the correctly reconstructed signal events, there is a possibility that the 
bachelor pion could be swapped with a pion from the af (1260) TT^TT^TT^ decay. We find 
that the probability of obtaining the wrong combination is negligibly small; the bachelor 
pion is unambiguously determined. 

We extract the signal yield by using d^E, Mbc and cos 6', where 6 is defined as the 
angle between the normal to the 05*^(1260) decay plane and the bachelor pion direction 
in the a^(1260) rest frame. This angle is employed to discriminate between the signal and 
B^ — > a^(1320)7r^ events, which have the same distributions in and Mbc as the signal 
events. 

For the correctly-reconstructed signal and B^ a^(1320)7r^ components, we use a sum 
of two bifurcated Gaussians and Gaussian functions to describe the AE and Mbc shapes, 
respectively. The possible differences between data and MC in the shapes of AE and Mbc 
distributions are taken into account using a — D~n^, D~ K^h'tt' decay control 
sample. The cos^ distributions are polynomial functions obtained from the MC samples. 
The Ai?-Mbc and cos 9 PDFs for SCF events are modeled by a smoothed two-dimensional 
histogram and a polynomial, respectively. 

The PDFs for charm 6 — >■ c and charmless h ^ u backgrounds are obtained from a 
large MC sample and modeled as two-dimensional smoothed histograms for Aii^-Mbc and 
polynomial functions for cos^. The charmless b —>■ u decay background contains a peaking 
background that has the same shape in A£^-Mbc as that of the signal. We estimate the 
peaking background yield using the MC simulation, and find that the dominant B decay 
modes are B^ K*~^7c~, Kq'^tt~ and For these modes 32 events are expected. The 



5 



uncertainty in this background component is included in the systematic error. 

We model the continuum background event shapes as a second-order polynomial and an 



ARGUS function for and Mbc, respectively. The parameters of the functions are 
determined from the fit. We model the cos^ distribution as a polynomial function. The 
parameters of the functions are obtained using events in the Mbc sideband region defined as 
Mbc < 5.26 GeV/c^. The parameters are fixed in the fit. 

We perform a 3D ( AE'-Mbc-cos ^) unbinned extended maximum likelihood fit to 31725 
candidates. The likelihood value is 

L = exp(-^n,)n[E^i^i(^^i'^bc,i,cos^i)], (1) 

where i runs over all the (1260)7r^ candidates, j indicates one of the following 

event categories: correctly reconstructed signal, SCF signal, the correctly reconstructed 
5° ^ a^(1320)7r^, SCF 5° a^(1320)7r^, b ^ c, b ^ u, and continuum backgrounds. 
rij is the yield of each category j, and Pj(A£', Mbc, cos^) is the PDF of the category j as 
a function of AE, Mbc and cos 9. The yields of the correctly reconstructed and SCF signal 
events are parameterized as (1 — /scF)^ai7r and /scF^^aivr, respectively, where /scf = 0.199 is 
the SCF fraction determined using the signal MC simulation, and Ua^^, is the signal yield. We 
fix the /scf value in the fit. The same parameterization is applied to the — > 0^(1320)7?"^ 
decay mode, where the SCF fraction (0.158) is also determined from the MC. In the 3D AE- 
Mbc-cos ^ fit, we vary the yields nj and parameters that determine the continuum background 
shape in AE and Mbc. The fit yields n^^^ = 654 ±70 signal events and n^aTr = 47±50 5° — 
(1320)7r''' events. Figure [1] shows the AE, Mbc and cos 6* distributions with projections of 
the fit results. 

We calculate the product of the branching fractions using 

^ af (1260)7r^)i3(af (1260) ^ n^n^n^) = ^ , (2) 

^BB ' ^dct ■ £^PID 

where e^et = 9.06% is the detection efficiency estimated using the signal MC simulation 
assuming a^(1260) p'^Tr^ decay via an 5- wave. The correction factor Epid = 0.90 takes 
into account the PID selection efficiency difference between the real data and the MC sim- 
ulation. It is obtained from a large sample of D*^ D^ir'^, K~n^ decays in data. 
We measure the product branching fraction to be (14.9 ± 1.6) x 10~^, where the error is 
statistical. 

To validate our results, we measure the branching fraction of the decay 5° D~7r~^, 
D- K+Ti'Ti- as a control sample. We obtain i3(5° D-7i+) = (2.42 ± 0.05) x 10-^ 
where the error is statistical only. The value is consistent with the World Average (W.A.) 
value of (2.54 ± 0.28) x 10"^ @. To estimate the fit bias, we perform fits to an ensemble 
of toy-MC pseudo-experiments and MC samples. The signal-to-background fractions used 
are obtained from the fit to the data. No bias is found in the fits to pseudo-experiment 
samples. We find a small difference (3.8 ±5.3%) for the fit to the MC samples, which is also 
consistent with zero. We assign a systematic uncertainty of ±3.8% to allow for this. 

We determine the af (1260) mass (M^J and width (P^J from an unbinned maximum 
likelihood fit to the af (1260) — * TT^TT^TT^ three pion mass distribution of the 2915 candidates 
in the signal box defined as \AE\ < 0.05 GeV and Mbc > 5.27 GeV/c^. We employ a 
relativistic Breit-Wigner (BW) function to model the (1260) mass distribution [l^ , which 
is multiplied by the mass dependent efficiency to form a PDF for the correctly reconstructed 
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FIG. 1: (a) projection in AE for events with Mbc > 5.27 GeV/c^, (b) projection in Mbc for 
events with \AE\ < 0.05 GeV, and (c) projection in cos 9 for events with Mbc > 5.27 GeV/c^ and 
\AE\ < 0.05 GeV. The points with error bars are the data. The soUd (dashed) curve shows the fit 
results, which includes the sum of all the components (backgrounds). 



signal events. We use a MC-determined PDF for the signal SCF events with the mass- 
dependent weight determined by the BW function. For the correctly reconstructed 5° — > 
a^(1320)7r^ events, we use a sum of three Gaussians. Since the B'^ a^(1320)7r^ SCF 
contribution is very small (< 0.1%) in the signal box, we neglect it. The PDFs of the charm 
and charmless B backgrounds are determined using the MC samples. The continuum PDF 
is obtained from the Mbc sideband events by subtracting the charm B decay contamination. 
In the fit, only two parameters, M^^ and r^^, are floated. The fractions of the signal to the 
backgrounds are fixed to the value determined from the 3D A£'-Mbc-cos^ fit. The fit yields 
Ma, = 1233^47 MeV/c^ and Ta, = 594:1l\11 MeV/c^ where the errors are statistical only. 
The measured af (1260) mass and width are comparable with the W.A. values [9i]. Figure [2] 
shows the three pion mass distribution along with the fit results. 

To estimate possible contributions from non-resonant B^ decays into p^vr+vr" and 
Tr+Tr'Tr+TT", we make use of the three pion mass distribution, since the 3D Ai?-Mbc-cos^ 
fit cannot distinguish between the signal and non-resonant components. We include PDFs 
for the non-resonant decay modes in the fitting likelihood function. The PDF of the non- 
resonant components is determined from large MC samples generated assuming three- and 
four-body phase space distributions; we use a single PDF for both modes since the shapes 
of these distributions are nearly the same. We divide the correctly reconstructed signal 
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FIG. 2: Three pion mass distribution with the fit result. The dashed and dotted curves show 
the signal and the sum of backgrounds, respectively. The solid curve indicates the sum of all the 
components. The points with error bars are the data. 

fraction in the signal box, /gig = 17.2%, into af (1260)7r^ and non-resonant contributions; 
the of (1260)7r^ and non-resonant fractions in the fit are parameterized as /sig(l — /nr) and 
/sig/nr, respectively, where /nr is the non-resonant fraction. The fraction /nr, and F^^ 
are allowed to float in the fit. We obtain /nr = 0.11±0.07, and M^^ and F^j values consistent 
with the results from the nominal fit where /nr = 0. 

The dominant systematic error of the branching fraction originates from the uncertainty in 
the non-resonant contributions, which corresponds to a ±11% systematic error. By varying 
the PDF parameters, we estimate the PDF systematic uncertainty of ±6.9%. We assign a 
±4.8% tracking reconstruction efficiency uncertainty. The systematic error on the detection 
efficiency mainly originates from the uncertainties in the af (1260) decays into p^vr^ via 
a D-wave. We assign a ±4.7% systematic error due to these uncertainties. We vary the 
branching fractions of the peaking backgrounds by one standard deviation, and add the 
differences from the nominal fit result in quadrature to determine systematic uncertainties 
of ±1.7%. Other sources of systematic error are the uncertainties in PID selection efficiency 
(±1.3%), fit bias (±3.8%) and number of BB pairs (±1.3%). We add each contribution in 
quadrature to obtain the total systematic error of ±15.3%. 

In summary, we measure the branching fraction product of B^ of (1260)%^ and 
05*^(1260) — >■ TT^TT^TT^ using a data sample containing 535 x 10® BB pairs: B{B^ 
af(1260)7rT)i3(af (1260) ^ Ti^n^n^) = (14.9 ± 1.6 ± 2.3) x 10-^ where the first and second 
errors are statistical and systematic, respectively. If we assume a 50% branching fraction for 

(1260) ^ 7r±7r±7rT decays, B{B'^ (1260)7rT) = (29.8 ± 3.2 ± 4.6) x 10"®. Our result 
is consistent with other measurements [5]. 
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